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The problem of anomalous diffusion in the momentum space is considered on 
the basis of the appropriate probability transition function (PTF). New general 
equation for description of the diffusion of heavy particles in the gas of the light 
particles is formulated on basis of the new approach similar to one in coordinate 
space [1]. The obtained results permit to describe the various situations when the 
probability transition function (PTF) has a long tail in the momentum space. The 
effective friction and diffusion coefficients are found. 
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I. INTRODUCTION 

Interest in anomalous diffusion is conditioned by a large variety of applications: semicon- 
ductors, polymers, some granular systems, plasmas in specific conditions, various objects in 
biological systems, physical-chemical systems, et cetera. 

The deviation from the linear in time < r 2 (t) >~ t dependence of the mean square 
displacement have been experimentally observed, in particular, under essentially non- 
equilibrium conditions or for some disordered systems. The average square separation of 
a pair of particles passively moving in a turbulent flow grows, according to Richardson's 
law, with the third power of time [2]. For diffusion typical for glasses and related com- 
plex systems [3] the observed time dependence is slower than linear. These two types of 
anomalous diffusion obviously are characterized as superdiffusion < r 2 (t) >~ t a (a > 1) 
and subdiffusion (a < 1) [4]. For a description of these two diffusion regimes a number 
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of effective models and methods have been suggested. The continuous time random walk 
(CTRW) model of Scher and Montroll [5], leading to strongly subdiffusion behavior, provides 
a basis for understanding photoconductivity in strongly disordered and glassy semiconduc- 
tors. The Levy- flight model [6], leading to super diffusion, describes various phenomena as 
self-diffusion in micelle systems [7], reaction and transport in polymer systems [8] and is 
applicable even to the stochastic description of financial market indices [9]. For both cases 
the so-called fractional differential equations in coordinate and time spaces are applied as 
an effective approach [10]. 

However, recently a more general approach has been suggested in [1], [11], which avoid 
the fractional differentiation, reproduce the results of the standard fractional differentiation 
method, when the last one is applicable, and permit to describe the more complicated cases 
of anomalous diffusion processes. In [12] these approach has been applied also to the diffusion 
in the time-dependent external field. 

In this paper the problem of anomalous diffusion in the momentum (velocity) space will 
be considered. In spite of formal similarity, diffusion in the momentum space is very different 
physically from the coordinate space diffusion. It is clear already because the momentum 
conservation, which take place in the momentum space has no analogy in the coordinate 
space. 

The anomalous diffusion in the velocity space is weakly investigated. Some attempts to 
investigate the influence of the long tails of correlation functions in velocity space have been 
done recently by W. Ebeling and M.Yu. Romanovsky (Contr. Plasma Physics, in print, 
private communication). The consequent way to describe the anomalous diffusion in the 
velocity space is, according to our knowledge, still absent. 

In Section II the diffusion equation in coordinate space for a homogeneous system is 
shortly reviewed. The diffusion in velocity space for the cases of normal and anomalous 
behavior of the probability transition function PTF is presented in Sections III, IV. 

II. DIFFUSION IN THE COORDINATE SPACE ON THE BASIS OF A 

MASTER-TYPE EQUATION 

Let us consider diffusion in coordinate space on the basis of the master equation, which 
describes the balance of grains coming in and out the point r at the moment t. The structure 
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of this equation is formally similar to the master equation in the momentum space (see, e.g., 
[1] , [11] ) . Of course, for coordinate space there is no conservation law, similar to that in 
momentum space: 

df g (r,t) 



dt 



J dv'{W(r,r')f g (T f ,t)-W(v',v)f g (r,t)}. (1) 



The probability transition W(r, r') describes the probability for a grain to transfer from 
the point r' to the point r per unit time. We can rewrite this equation in the coordinates 
p = r' — r and r as: 

df g (r,t) 



dt 



j dp {W(p, r + p)f g (r + p,t)- W(p, r)f g (r, t)} . (2) 



Assuming that the characteristic displacements are small one may expand Eq. (j2J) and arrive 
at the Fokker-Planck form of the equation for the density distribution f g (r,t) 



df g (T,t) _ d 



dt dr r 



A a (r)f g (r,t) + l-(B aP (r)f g (v,t)) 



(3) 



The coefficients A a and B a p, describing the acting force and diffusion, respectively, can be 
written as functionals of the PTF in the coordinate space W (with the dimension s) in the 
form: 

A a (r) = J d s p Pa W(p,r) (4) 

and 

B a p{r 



~ j d s pp aP pW(p,r). (5) 



For the isotropic case the probability function depends on r and the modulus of p. For a 
homogeneous medium, when r-dependence of the PT is absent, the coefficients A a = while 
the diffusion coefficient is constant with B a p = 5 a pB, where B is the integral 

B = hS dSpp2w ^- (6) 

This consideration cannot be applied to specific situations in which the integral in Eq. ([6]) 
is infinite. In that case we have to examine the general transport equation (CQ). We will now 
consider the problem for the homogeneous and isotropic case, when the PT function depends 
only on \p\. By Fourier-transformation we arrive at the following form [11] of Eq. ([T]): 

df g (k,t) 



dt 



j <r P [exp(ikp) - l]W(\p\)f g (k,t) =X(k)/ fl (k,t), (7) 
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where X(k) = X(k). Let us assume a simple form of the PT function with a power depen- 
dence on the distance W(p) = C/\p\ a , where C is a constant and a > 0. Such type singular 
dependence is typical for jump diffusion probability in heteropolymers in solution (see, e.g., 
[13], where the different applications of anomalous diffusion are considered on the basis of 
the fractional differentiation method). For the one- dimensional case we find: 

X(k) = -A^dustn 2 f^j W{u) = -2 3 - a C\k\ a ~ l jT ' d(^^. (8) 

For the values 1 < a < 3 this function is finite and equal to 

Cr[(3- a)/2]\k\^ 
X{k) ~ -2%^r(a/2)(a-l)' { ) 

where T is the Gamma-function. At the same time the integral in Eq. ([6]) for such a type of 

PT functions is infinite, because usual diffusion is absent. 

At the same time the integral ([6]) for such a type of PT functions is infinite, because usual 

diffusion is absent. The procedure considered for the simplest cases of power dependence of 

the PT function is equivalent to the equation with fractional space differentiation [10], [13]: 

^^ = CA^f g (x,t), (10) 

where A M//2 is a fractional Laplacian, a linear operator, whose action on the function f(x) in 
Fourier space is described by A^/ 2 /(x) = — {k 2 )^ 2 f(k) = — \k\ tM f(k). In the case considered 
above p = (a — 1), where < p < 2. For more general PT functions, which (for arbitrary 
values p) are not proportional to the a power of p, the method described above is also 
applicable, although the fractional derivative does not exist. 

For the case of purely power dependence of PT the non-stationary solution for the density 
distribution describes so-called super-diffusion (or Levy flights). The solution of Eq. ffTUl) in 
Fourier space reads: 

f g (k,t)=exp{-C\k\' i t), (11) 
which in coordinate space corresponds to a so-called symmetric Levy stable distribution: 

1 , 



For the general case it follows from Eq. (|7j) that 



x n 
/i,0 



(12) 



f g (k,t) = C 1 exp[X(k)t], (13) 
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with some constant C\. 

The consideration on the basis of PTF function given above, permits to avoid the frac- 
tional differentiation method and to consider more general physical situations of the non- 
power probability transitions for arbitrary space dimension. 



III. DIFFUSION IN THE VELOCITY SPACE ON THE BASIS OF A 

MASTER- TYPE EQUATION 

Let us consider now the main problem formulated in the introduction, namely, diffusion 
in velocity space (V-space) on the basis of the respective master equation, which describes 
the balance of grains coming in and out the point p at the moment t. The structure of this 
equation is formally similar to the master equation in the coordinate space Eq. (j2j) 

^f^ = J dq{W(q l p + q)/ a (p + q 3 t)-W(q l p)/ p (p l t)}. (14) 

Of course, for coordinate space there is no conservation law, similar to that in the momen- 
tum space. The probability transition lF(p, p') describes the probability for a grain with 
momentum p' (point p') to transfer from this point p' to the point p per unit time. The 
momentum transferring is equal q = p' — p. Assuming in the beginning that the character- 
istic displacements are small one may expand Eq. (j2J) and arrive at the Fokker-Planck form 
of the equation for the density distribution f g (p,t) 



df g (p,t) d 



(15) 



dt dp, 

A a {p) = J d s qq a W(q,p); B a(3 (p) = d s qq a q p W(q,p). (16) 

The coefficients A a and B a p describing the friction force and diffusion, respectively. 

Because the velocity of heavy particles is small, the p-dependence of the PTF can be 
neglected for calculation of the diffusion, which in this case is constant B a p = 5 a pB, where 
B is the integral 

B = Ysj dS ^ w ^- 

If to neglect the p-dependence of the PTF at all we arrive to the coefficient A a = 
(while the diffusion coefficient is constant). This neglecting, as well known is wrong, and the 
coefficient A a for the Fokker-Planck equation can be determined by use the argument that 
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the stationary distribution function is Maxwellian. On this way we arrive to the standard 
form of the coefficient MTA a (p) = p a B, which is one of the forms of Einstein relation. For 
the systems far from equilibrium this argument is not acceptable. 

To find the coefficients in the kinetic equation, which are applicable also to slowly de- 
creasing PT functions, let us use a more general way, based on the difference of the velocities 
of the light and heavy particles. For calculation of the function A a we have take into account 
that the function W(q, p) is scalar and depends on q, q • p,p. Expanding W(q, p) on q ■ p 
one arrive to the approximate representation of the functions W(q, p) and W(q, p + q): 

W(q, p) ~ W(q) + W'(q){ci • p) + ±W"(q)(q ■ p) 2 . (18) 

W(q, p + q) ~ W{q) + W\q) (q ■ p) + ^W"{q)(q ■ p) 2 + q 2 W'(q), (19) 

where W'{q) = dW(q, q • p)/9(qp) | q . p = and W"{q) = d 2 W(q, q • p)/<9(qp) 2 | q . p=0 . 
Then, with the necessary accuracy, A a equals 

A a {p) = J d s qq a q m W'(q)=p a J d s qq a q a W'{q) = ^- J d s qq 2 W'(q) (20) 

If for the function W(q,p) the equality W'(q) = W(q) jlMT is fulfilled, then we arrive to 
the usual Einstein relation 

MTA a (p)=p a B (21) 

Let us check this relation for the Boltzmann collisions, which are described by the PT- 
function W(q, p) = w B (q, p) [11]: 



2n [°° da 
H Q Jq/2f, do 



arccos (1 - - - ),u 
2fi z u z 



f b (u 2 + v 2 - q -v/fx), (22) 



where (p = Mv) and da /do and /& are respectively the differential cross-section for scat- 
tering and the distribution function for the light particles. For the equilibrium Maxvellian 
distribution /° the equality W'(q) = W(q)/2MT is evident and we arrive to the usual 
Fokker-Planck equation in velocity space with the constant diffusion D = B/M 2 and fric- 
tion {3 = B/MT = DM/T coefficients, which satisfy the Einstein relation. 

For some non-equilibrium situations the PTF can possess a long tail. In this case we have 
derive a generalization of the Fokker-Planck equation in spirit of the above consideration for 
the coordinate case, because the diffusion and friction coefficients in the form Eqs. ffTTl) . (|20l 



7 



diverge for large q if the functions have the asymptotic behavior W(q) ~ l/q a with a < s + 2 
and (or) W'(q) ~ 1/q* 3 with (3 < s + 2. 

Let us insert in Eq. (fl4l) the expansions for (as an example we choose s = 3, the 
arbitrary s can be considered by the similar way). With necessary accuracy we find 

dfg(p,t) _ 



dt 



dq{f g (p + q, t) [1 + q a d/dp a ] [W(q) + W'(q) (q • p) + 



\w"(q)(q ■ p) 2 ] - f g (p, t)[W(q) + W\q) (q • p) + ^"(g)(q • p) 2 ]} (23) 
After the Fourier-transformation /(r) = / ^^-exp(ipr)/(p, i) Eq. (123]) reads: 



/ dq{exp(-i(qr) [W(g) - iW (g) (q • — ) + 



~W"{q){q ■ ^f) - [W(q) - iW'(q) (q ~) - ^"(g)(q ■ |;) 2 ]}/ 3 (r, f) (24) 
We can rewrite this equation as 



d 4M = A(r)f{r) + fl.(r) A /(r , t) + C£ ,, (r) _|_ /(r , ( ) (25) 



where 

/* r°° r m i /it i 

^(?) (26) 



//"OC 
dq[exp(-z(qr)) - l]W(g) = Air dqq 2 



POO 

/ dgg 2 


sin (qr) 


/o 


qr 



B a = r a B(r); B(r) 



— / dqqr[exp(— z(qr) — l]W"'(g) 



47T 



^ 2 -/o 



dgg 2 



cos (qr) — 



sin(qr) 
qr 



W'{q) 



(27) 



C a p{r) = r a r p C(r) = -- / dqg a g^[exp(-i(qr) - l]W"(g) 



(28) 



2tt 



C(r) 



cigg 4 



1 

2H 



dq(qr) 2 [exp(— i(qr) — l]iy"(g) 



2sin(qr) 2cos (qr) sin(qr) 1 

qr 3 



W"(g) 



(29) 



For the isotropic function /(r) = /(r) one can rewrite Eq. (I25I) in the form 



df g (r,t) 
dt 



d 



d 2 



A(r)f(r) + B(r)rj-f(r) + C(r)r 2 j-^f(r) 



(30) 
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For the case of strongly decreasing PDF the exponent under the integrals for the functions 
A(r), B(r) and C(r) can be expanded 



Air) ~= -— / dqq 2 W(q); B(r) 
6 " 



-- / dqq 2 W'(q); C(r) ~ 0. 
6 



(31) 



Then the simplified kinetic equation for the case of short-range on q- variable PTF (non- 
equilibrium, in general case) reads 



d ^j^ = A r 2 f(r)+B r^-f(r) 

where A = -1/6 / dqq 2 W(q) and B = -1/6 / dq,q 2 W'(q). 
Stationary solution of Eq. ( 1301) for C{r) =0 reads 

j Mr') 



(32) 



f g (r,t) = Cexp 



dr' 



r'Bir' 



Cexp 



Apr' 
2B n 



The respective normalized stationary momentum distribution equals 

B p\ 



fM 



N R 3//2 



(2nA y/i eXP ^ 2A 

Therefore in Eq. (I3T|) the constant C = N g . Equation (1321) and this distribution are the 
generalization of the Fokker-Planck case for normal diffusion on non-equilibrium situation, 
when the prescribed W(q, p) is determined, e.g., by some non-Maxwellian distribution of 
the small particles /&. To show this by other way let us make the Fourier transformation of 
ff25|) with C = and the respective A and B a : 

dfg(p,t) 



(33) 



(34) 



dt 



d 2 d 
-A —f 9 (p,t)-B — Pa f g (p,t), 



(35) 



Therefore we arrive to the Fokker-Planck type equation with the friction coefficient (3 = —Bq 
and diffusion coefficient D = —Aq/M 2 . In general these coefficients (Eq. (I3T]) ) do not satisfy 
to the Einstein relation. 

In the case of equilibrium iy-function (e.g., /& = /°, see above) the equality W'(q) = 
W(q)/2MT b is fulfilled. Then A{r)/rB(r) = MT b r (A = MT b B ). Only in this case 
the Einstein relation between the diffusion and friction coefficients exists and the standard 
Fokker-Planck equation is valid. 



IV. THE MODEL OF ANOMALOUS DIFFUSION IN V - SPACE 



Now we can calculate the coefficients for the models of anomalous diffusion. 
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In this paper we calculate only the simple model system of the hard spheres with the 
different masses m and M 3> m, da /do = a 2 /A. Let us suppose that in the model under 
consideration the small particles are described by the prescribed stationary distribution 
fb = rib(fib/uo (where 4>b is non-dimensional distribution, uq is the characteristic velocity for 
the distribution of the small particles) and £ = (u 2 + v 2 — q ■ v//x)/uq. 

W a (q,p) = £p-r dC-<f> b (0- (36) 

If the distribution </>&(£) = l/£ 7 (7 > 1) possess a long-tail we get 



W (a v) - ^ ^ 7 r - Uba2n & (37) 
2// 2 u g (1 - 7) 2/i 2 u g (7 - 1) 



where £ = (? 2 /4/x 2 + w 2 — q ■ v/ h)/uq. 

For the case p = the value £ ~~ * £0 = 9 2 /4/i 2 Wg an d we arrive to the expression for 
anomalous = W a 

W a (q, p = 0) = ^ „ 2 = 4^7- (38) 

a ^ ' 2 3 -^( 7 -l)/i 4 - 2 ^- 2 V 7 - 1 

The function A(r), according to Eq. (1261) 

sin(qr) 



A(r) 



/■oo 

47r / rfgg 2 


sin (qr) 


Jo 


qr 



poo 

W(q) = 4irC a / dq 
Jo 



(39) 



qr 

Comparing the reduced equation (see below) in the velocity space with the diffusion in 
coordinate space (27 — 1 <-» a and W(q) = C/q 21 " 1 ) we can establish that the convergence 
of the integral in the right side of Eq. fl39|) (3d case) is provided if 3 < 27 — 1 < 5 or 
2 < 7 < 3. The inequality 7 < 3 provides the convergence for small q (q — > 0) and the 
inequality 7 > 2 provides the convergence for q — ► 00. 

Now to determine the structure of the transport process and the kinetic equation in the 
velocity space we have find the functions W'(q) and W"(q). 

If p ^ to find W'(q) and W"(q) we have use the full value £ = (<? 2 /4/i 2 +p 2 /M 2 - 
q ■ p/M(/,)/uq and it derivatives on q • p at p = 0, £q = — 1/M//W 2 , and £q = 0- Then 

\xrn \ nba2lT t-7 u/"/- ^ t - 7 -i ,„ n x 

^ (q ' p) "2M^^ o; ^ (q ' p) "2MV^ e ° (40) 

Therefore for p = (£0 - ► £0) we obtain the functions 

W'(a)- (V^) 7 ^ . T ^ fq) - (V^) 7+ 
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We have establish now the conditions of convergence the integrals for B(r) and C(r). 



B(r) 



% I dqq 2 



cos (qr) 



sin(qr) 
qr 



W'(q) 



(42) 



Convergence B(r) exists for small q if 7 < 2 and for large q —>■ 00 for 7 > 1/2. 
Finally for C(r) convergence is determined by the equalities 7 < 2 for small q and 7 > 1 
for large q 



CM = % 



dqq 4 



2sin(qr) 2cos (qr) sin(qr) 1 



W"(q) 



(43) 



Therefore to provide convergence for A, B, C for large q we have provide convergence for 
A, that means 7 > 2. To provide convergence for small q enough to provide convergence for 
B and C, that means 7 < 2. Therefore for the purely power behavior of the function /&(£) 
convergence is absent. However, for existence of the anomalous diffusion in the momentum 
space in reality the convergence for small q is always provided, e.g. by finite value of v or by 
change of the small g-behavior of W(q) (compare with the examples of anomalous diffusion 
in coordinate space [1]). Therefore the "anomalous diffusion in velocity space" for the power 
behavior of W(q), W'(q) and W"(q) on large q exists if for large q the asymptotic behavior 
of W(q — ► 00) ~ 1/q 2 " 1 ^ 1 with 7 > 2. At the same time the expansion of the exponential 
function in Eqs. (|26|) - (|29|) under the integrals, which leads to the Fokker-Planck type kinetic 
equation is invalid for the power-type kernels W{<\, p). 



V. CONCLUSIONS 



In the previous sections we shortly reviewed the anomalous diffusion in the coordinate 
space and firstly consequently considered the problem of anomalous diffusion in momentum 
(velocity) space. The new kinetic equation for anomalous diffusion in velocity space is 
established. For the normal diffusion the friction and diffusion coefficient are found for the 
non-equilibrium case. For equilibrium case the usual Fokker-Plance equation is reproduced 
as the particular case. The model of anomalous diffusion in velocity space is described on 
the basis of the respective expansion of the kernel in master equation and the conditions of 
the convergence for the coefficients of the kinetic equation are found. 
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